Density-functional theory of quantum wires and dots in a strong 

magnetic field 



T.H. Stoof and Gerrit E.W. Bauer 
Department of Applied Physics and Delft Institute for Microelectronics and Suhmicrontechnology, 

• Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands 

: (February 1, 2008) 

C ■ 

^3 ; 
>—i ; 

00 ■ Abstract 
^ ! 

, We study the competition between the exchange and the direct Coulomb 

m : 

, interaction near the edge of a two-dimensional electron gas in a strong mag- 

\ netic field using density-functional theory in a local approximation for the 

\ exchange-energy functional. Exchange is shown to play a significant role in 

^ \ reducing the spatial extent of the compressible edge channel regions obtained 

' from an electrostatic description. The transition from the incompressible 

Q \ edge channels of the Hartree-Fock picture to the broad, compressible strips 



X 



predicted by electrostatics occurs within a narrow and experimentally acces- 



H ' sible range of confinement strengths. 



PACS numbers: 73.20.Dx; Tl.lO.+x, 71.45. Gm 



Submitted to Phys. Rev. B. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The concept of current- carrying edge channels accounts for the magnetotransport prop- 
erties of a two-dimensional electron gas (2DEG) in a high magnetic field, both in the integeilll 
and fractional!^ quantum Hall regime. Although the initial theoretical studies have used a 
noninteracting picture of edge channels, a considerable effort has recently been devoted to un- 
derstand the effects of electron-electron interactions in the integer quantum Hall regime.il"i 

At present there are two incompatible pictures for the electronic ground state of edge 
channels in the integer quantum Hall regime. In the Hartree-Fock approximation the ground 
state wave function is a single Slater determinant, which corresponds to occupation numbers 
zero and one.i In this description the edge state is incompressible and the electron density 
drops rapidly at the edges, on a length scale which is typically of the order of the magnetic 
length. Although the Hartree-Fock approximation is widely used for the description of 
quantum dots and wires,illi it has been challenged on the grounds that is does not take the 
global electrostatics into account properly.Hll In the electrostatic description given in Ref. |^ 
a more gradual variation of the electron density at the edge is found energetically favorable. 
Here the typical length scale for the density drop at the edge is of the order of the depletion 
length, which is much larger than the magnetic length. In this picture the electrons at the 
edge can screen the electrostatic confinement potential by a slow density variation and the 
edge states are therefore compressible. However, this description completely neglects the 
exchange interaction. 

A transition from incompressible to compressible edge states was already qualitatively 
discussed in Ref. ^, where a spontaneous transition from an unpolarized {i.e. equal oc- 
cupation for both spin levels) to a polarized (different occupation for different spin levels) 
Hartree-Fock ground state was found for a critical confinement strength. It was speculated 
that a development to the electrostatic regime takes place as the confinement strength is 
decreased. More recently the transition between smoothly and abruptly varying density 
distributions has been studied by Chamon and Wen.i On the basis of few particle exact 
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solutions they predict formation of compressible edge states when the strength of the con- 
finement potential is reduced beyond a certain point. Their approach is difficult to extend 
to quantum wires or dots with a large number of electrons. 

In the present work we interpolate between the different regimes of applicability of both 
Hartree-Fock and electrostatic pictures by extending the Thomas-Fermi approach used in 
Ref. P to a Tho mas- Fermi- Z?2rac like treatment of the exchange effects in strong magnetic 
fields. First results of the present approach have been published in Ref. Very recently 
Ferconi and Vigna have studied the ground state energies and densities of a quantum dot 
in an arbitrary magnetic field, taking into account exchange-correlation effects by employing 
a Kohn-Sham scheme of current-density functional theory.lll For small quantum dots (2 or 
3 electrons) their method yields an accuracy better than 3% when compared with exact 
results. For a larger number of electrons and high magnetic fields we find that the Kohn- 
Sham formalism encounters serious problems. 

In Sec. |I| we describe our implementation of density-functional theory in strong magnetic 
fields as applied to quantum wires and derive expressions for the density profile and single- 
particle potentials of the ground state. In Sec. |T| we present results of the numerical 
calculations and we investigate both the accuracy of our theory and the relation with possible 



experiments. Two experimentally relevant applications are studied in Sec. jl^; the infiuence 
of a plane of constant potential parallel to the 2DEG on the electronic ground state and the 
ground state properties of a quantum dot in a strong magnetic field. Section summarizes 
our conclusions. 

II. DENSITY-FUNCTIONAL THEORY IN STRONG MAGNETIC FIELDS 

A. The system 

Let us consider first a quantum wire of the strictly two dimensional electron gas along 
the y axis in the xy plane (Fig. [l|). Perpendicular to the plane a strong uniform magnetic 
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field is applied. Also present is a uniform positive background charge which ensures global 
neutrality. We use periodic boundary conditions in the y-direction and adopt the Landau 
gauge so that Aq = Bqxj and Bq = B^z. In the high magnetic field limit considered here 
only the lowest (spin-up and spin-down) Landau levels are occupied. 

The wire is confined in the x direction by a parabolic confinement potential: 

v,{x) = « ^ x^ (1) 

where e = eoCj. is the dielectric constant and a is a dimensionless parameter. We will ignore 



here the difference between the 'bare' magnetic length Ib = yh/eBo and the magnetic length 
/ which is renormalized by the parabolic confinement potential: 



(2) 



where ujc = eBo/m* is the cyclotron frequency with m* the effective electron mass and 
cuq characterizes the parabolic confinement potential ^m*LjQx'^. Since we treat the problem 
numerically, the assumption of a parabolic form for the confinement is not essential: It is 
chosen here because it is widely used in the literature. All lengths will be given in units of 
Ib throughout the paper unless otherwise indicated. 



B. Theory 

We use density-functional theory in strong magnetic fields to find the electronic ground 
state of the system.0 Without a magnetic field, the total energy of the system would be a 
functional of the electron density only. However, a magnetic field, apart from giving rise to 
Zeeman splitting, causes orbital currents to flow in the electron gas, even when the system 
is in thermodynamic equilibrium. As a consequence, the total energy is now a functional of 
the density r;,(r), the spin density s(r) and the paramagnetic current density jp(r).lli 

We do not take into account correlation effects, which means that we operate strictly in 
the integer quantum Hall regime, since correlation is responsible for the energy gaps that 
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cause the fractional quantization. Due to the fact that the correlation part of the energy is 
disregarded and exchange does not depend on the current density the total energy functional 
depends only on the total and the spin density or equivalently, since the spin-quantization 
axis is well defined by the strong magnetic field, on the density of spin up and spin down 
electrons. 

In the Kohn-Sham scheme of density-functional theory the ground state densities n'^ of 
the interacting electron system are expressed in terms of a set of M'^ Kohn-Sham orbitals 
0f(r): 

^'^(r) = ElCWr, (3) 

i=l 

where cr =| , | denotes up or down spin respectively and M" is the number of electrons with 
spin a. These orbitals satisfy the Kohn-Sham equations: 

{^V^ + \/w}c(r) = 6^C(r), (4) 

with V"g'|j(r) the effective one-particle potential: 

V^^{v) = K(r) + VniHr) + i[n\ n^]; r) , (5) 

where Vh is the Hartree potential and the exchange potential. Here we have disregarded 
the exchange-correlation vector potential A^c which gives only a very small contribution to 
the total energy.0 The Kohn-Sham ground state wave function is a Slater determinant of 
the M°" lowest Kohn-Sham orbitals: 

$'^(ri, r2, VM^) = y^det {0J (r,)} . (6) 

It can be shown (see Sec. [IIIBD that in the limit of high magnetic fields the mixing 



between Landau level wave functions vanishes. Consider, e.g., the case of a single occupied, 
spin resolved Landau level. For large magnetic fields the Kohn-Sham ground state wave 
function, Eq. (^), reduces to a Slater determinant of lowest Landau level wave functions 
with integer filling. As a consequence the scheme cannot describe fractional filling in the 



extreme quantum limit we are interested in. The reason for the failure of this procedure is 
the fact that the crucial assumption of the Kohn-Sham scheme, namely that the interacting 
f -represent able densities are also non-interacting f-representable,Ei breaks down in the high 
magnetic field limit. 

It is not possible to extend the Kohn-Sham scheme to fractional filling since in that case 
the effective potential, ^^'^^(r), is no longer a unique functional of the total and spin density. 
However, at sufficiently high magnetic fields the basis wave functions are known to be just 
Landau level wave functions. We may then use a variational scheme which allows fractional 
filling and in which E[n\n^] is a unique functional of the total and spin density. 

The equilibrium density distribution n{x) = n\x) + n^{x) which minimizes the total 
energy E of the system can be found by solving: 

6E[n\n^] _ 6E[n\n^] _ 
Sn^ix) ~ Sn^ix) ~ ^' ^ ^ 

where n is the chemical potential and the total energy of the quantum wire is given by: 

E[n\n^] = T[n] + E^[n] + E^[n] + Ez[n\n^] + E^[n\n^]. (8) 

The total energy, Eq. (|^), consists of five contributions: The first term on the right hand 
side is the kinetic energy which we define as a functional of the total density: 

/oo 
dxn%ix) iN + -)hu;,, (9) 

where labels the Landau levels. The densities are the partial densities for a given 

spin direction a and Landau level N. The total electron density is found by summing over 
all occupied Landau levels and over spin directions: n{x) = J^Na ^Ar(2^) and is for sufficiently 
high magnetic fields given by: 

^i^) = J2 '^Nx \^Nx{x,y)\^, (10) 

uNX 

where X = kyl\ is the quantum number of an electron with momentum ky, ipNx{x,y) ~ 
4>Nx{x) exp {iXy/l\) are the single particle bulk Landau level wave functions and v'}qx is the 
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local filling factor for electrons with spin a in Landau level A^. For well behaved confinement 
potentials the partial densities can be deduced from the total densities and For these 
potentials the total filling factor for a given spin direction, z/'^(x), lying between Nmax — 1 
and Nmax, always consists of N^ax — 1 completely filled Landau levels and a partially filled 
one. This means that the kinetic energy (and also the exchange energy, see Eqs. (|T3p and 
(|T6|)) is still a functional of the total and not the partial densities. 

The second contribution is the confinement energy which is given by: 



Ec[n] = / dx n{x) eVc{x). (11) 

J —oo 

The electrostatic Hartree energy is: 

/oo 
dx n{x) eVnin], (12) 
-oo 

where the Hartree potential, Vh [n] , is itself a functional of the total density. For the strictly 



two dimensional electron gas (See Ref. |T4| for the case of a quasi-two dimensional electron 
gas) it is given by: 



Vff[r;,] = — / dx n{x') ln|a; — a;'|. (13) 

2nelB J-oo 

The fourth term denotes the Zeeman energy: 

Ezin' ,n^] = -q^ibBq / dx In' {x) - {x) \ , (14) 

Z J — oo ^ ^ 

where g is the bare Lande factor. 

Our only concern left is the explicit form of the last term in Eq. (^, which is the 
exchange-energy functional. In the local density approximation (LDA) it reads: 

Ex[n\n'^] = / dx \'n){x)ex{n'{x))+n'^{x)ex{n^{x))\. (15) 

J —oo ^ ^ 

Here denotes the one-particle exchange energy of the homogeneous 2DEG with ground 
state densities and n^, which is magnetic field dependent: In the absence of a magnetic 
field, it is proportional to y/n (Ref. |T3p, but in the magnetic quantum limit considered here, 
the local exchange energy per electron is: 



Gxin^ix)) = -- — — J2 ^NN' n%,{x), (16) 

where the coefficients c^n' describe the exchange couphng between Landau levels N and 
A^'. They can be found by calculating the exchange energy per particle of an extended 
Landau level system.lli For two occupied levels, e.g., they are cqo = V^tt^, cqi = cio = y^7r'^/2 
and cii = y^OTT'^/S. 

The local density approximation is clearly only justified when the density variations are 
small on a characteristic length scale, which in our case is the magnetic length. In practice, 
however, it often turns out to be a useful and accurate tool even in cases in which this 
condition is not met. El 

Because we expect the results for higher Landau levels to be qualitatively the same, we 
will in the following restrict ourselves to the magnetic quantum limit, for which only the 
lowest Landau level is occupied. In the space of the lowest (spin-polarized) Landau level, 
the kinetic energy per particle is constant and may be disregarded. The validity of this 



approximation is discussed in Sec. |III B| . The density is in this case given by: 



n{x) = -^E r z.'^(X) e-(^-^)\ (17) 

where we have replaced the sum over X in Eq. (|l^) by an integral and substituted the lowest 
Landau level wave functions. This relation shows that in the extreme quantum limit there 
exists a one to one correspondence between the density and the ffiling factor which enables 
us to use the latter as the variational function in the determination of the ground state of 
the system (see Sec. pi A| ). 



The exchange energy per electron in the lowest Landau level (A^ = A^' = 0) reduces to: 



i.e. the one-particle exchange energy for a given spin direction is just proportional to the 
corresponding electron density. 

On substituting Eqs. (|T3|), ( P!7| ) and (|TB|) in Eq. (g) for the total energy and performing 
the functional derivatives, Eqs. (|^) become: 
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8E[n\n^] 
5E[n\n^] 



eVc{x) + eVnin] + 2e^[n^] - l-gfiBBo = ii. 



(19b) 



(19a) 



The numerical treatment of these coupled equations and the results will be presented in the 
next section. 

III. CALCULATION OF THE GROUND STATE PROPERTIES 

A. Results 

As described above we obtain the ground state occupation numbers by numerically min- 
imizing the total energy with respect to the filling factors, while keeping the total number 
of electrons constant. We discretize the filling factor and impose the boundary conditions: 
< z/'^ < 1. The numerical algorithm that minimizes the total energy uses a sequential 
quadratic programming method and is very stable. Numerical integration of the singular 
integrand in Eq. ( [TB| ) requires some care. By sampling the integrand equidistantly and sym- 
metrically around the divergence, x = x', and by sampling enough points to avoid oscillatory 
behavior in the resulting integral. An advantage of the present numerical scheme is that it 
can be very easy generalized to different confinement potentials. It should be noted, how- 
ever, that the choice of confinement potential is not completely arbitrary. For e.g. a hard 
wall potential there is a stronger mixing with higher Landau levels near the edges and for 
such a confinement Eqs. {^i7\ ) and (0) may be no longer accurate. 

With the approach outlined above we can obtain, to a very good approximation, both 
the electrostatic solution (by neglecting the exchange term) and the Hartree-Fock solution as 
calculated by Dempsey et a/.i (by forcing the filling factors to be integer valued). Therefore 
we expect our approach to give a good description of the intermediate regime between these 
two extremes. 

In Fig. 1^ ground state occupation numbers obtained with different methods are plotted 
as a function of the confinement strength a. Fig. ^(a) shows the purely electrostatic solution. 
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the Hartree-Fock solution is plotted in Fig. ^(b) and the solution including exchange in the 
local density approximation is shown in Fig. @(c). The small incompressible regions found 
in the electrostatic solution are caused by Zeeman splitting. 

The calculations were performed using a magnetic field of 7.2 T which for filling factor v = 
2 in the bulk corresponds to a zero-field bulk density of 3.5x10^^ cm^^. For GaAs the static 
dielectric constant = 12.5 and the bare Lande factor \g\ = 0.44. The filling factors are 
plotted for a constant number of electrons and confinement strengths a = 0.035, 0.041 and 
0.047 respectively. Converted to energy level spacings of a parabolic confinement potential 
of the form: ^m*ujQx'^ these values correspond to Hujq = 3.4 meV, 3.8 meV and 4.1 meV 
respectively. 

As can be seen from Fig. |](c) the effect of exchange on the electrostatic solution is a 
reduction of the compressible regions because it favors integer filling. However, a comparison 
with Fig. H(b) shows that for soft confinement potentials the Hartree-Fock approximation no 
longer holds, i.e. the electrostatic interaction overcomes the tendency of exchange to form 
an incompressible ground state. Note, however, that the solution for a = 0.047 in Fig. ^(c) 
is almost identical to the corresponding Hartree-Fock solution shown in Fig. 0(b), which 
shows that for confinement potentials that are strong enough the solution is forced into the 
integer filling regime, in agreement with the qualitative picture given in Ref. |. Note that 
this transition from incompressible to compressible state is a genuine correlation effect since 
it corresponds to the mixing of different single-determinant configurations, in spite of our 
exchange-only potential. 

It is clear from Fig. ^(c) that the width of the incompressible region is of the same order 
of magnitude as the width of the outermost compressible strip, even for the relatively wide 
strips we consider here (approximately 420 nm for a magnetic field of 7.2 T). This in contrast 
to the electrostatic description of edge states shown in Fig. 0(a), where the width of the 
incompressible region is always much smaller than that of the compressible one. However, 
for soft confinement potentials the solution including exchange does not deviate dramatically 
from the purely electrostatic one, which obviously implies that the Hartree-Fock approach 
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breaks down in this regime. Only a strong confinement potential can cause incompressible 
edge channels, although in Fig. ^(c) the confinement is still not strong enough to reduce 
the splitting to that of the bare Zeeman splitting. In the case of even harder confinement 
[a ~ 0.073, i.e. fiujo ~ 5.3 meV) we find that the splitting indeed reduces to this minimum 
value (not shown). 

Note that the entire range from the electrostatic to the (unpolarized) Hartree-Fock regime 
can be realized in realistic confinement potentials with Uujq ranging from 3.4 meV to 5.3 meV. 
Confinement potentials with level spacings of a few meV have been realized0lli so it should 
be possible to test these results by experiments. A possible method to stiffen the confine- 
ment experimentally is fabricating the 2DEG closer to the surface of the AlGaAs/GaAs 
heterostructure. 

In essence, our method to include exchange is using an position dependent (^-factor that 
depends on the local density. A simple approximation to this method would be to use the 
enhanced gf-factor corresponding to a single occupied spin leveL^ This method does increase 
the width of the incompressible strip but in general does not reproduce the qualitative and 
quantitative features obtained with our method.0 

In Fig. ^ the calculated one-particle potentials in the wire are plotted for both spin 
directions. These potentials correspond to the ground state density distribution of Fig. ||(c) 
with a = 0.035 which is also included in the figure. Fig. |^(a) shows the electrostatic 
potential, which consists of the Hartree and confinement potential. Also plotted are the 
exchange potential and the total potential for the majority-spin electrons, i.e. electrons 
occupying the lowest spin level. In Fig. ^(b) the same can be seen for the minority-spin 
electrons. In the figures the constant Zeeman term has been omitted for simplicity so that 
the total potential is just the sum of the electrostatic and exchange potential. 

Fig. ^(a) shows that for a majority spin electron the energetically most favourable po- 
sition is near the edges. This is because these electrons are forced by exchange to form 
an electrostatically unfavourable density profile in which no screening of the confinement 
potential is possible, i.e. the density of majority spin electrons is constant in the bulk of 

11 



the wire. This in contrast to the minority spin electrons, for which the potential in the bulk 
is fiat. This is due to the fact that these electrons still have the freedom to form a density 
distribution which can screen the external potential without being influenced much by the 
exchange interaction. 

The typical voltage drop from the edge to the middle of the wire has experimentally 
realistic values of ~ 300 mV for the bare confinement potential and ~ 30 mV for the 
selfconsistent potential. 

We also performed calculations of the widths of compressible and incompressible regions 
as a function of the confinement strength. The results are visible in Fig. ^ where the widths 
of the outermost incompressible region (I) and innermost compressible region (C) have been 
plotted as a function of a. If no incompressible strip is present in the middle of the wire, C 
is defined as half the width of the total compressible middle region (see inset). In the figure 
the points do not always lie on the smooth curve which serves as a guide for the eye. This is 
due to the spacing between sampling points (0.2 Z^), which imposes a upper bound on the 
accuracy of the calculated widths. 

The compressible strip shrinks rapidly as the confinement is increased but if the width has 
reached a value of roughly a few magnetic lengths, the further decrease becomes very slow 
and in practice we always find a small but finite compressible region. It is clear, however, 
that a Hartree-Fock treatment should give good results for hard confinements where small 
compressible regions exist. We recover the spin-polarizing transition of Ref. | since the 
incompressible strip between spin-up and spin-down channels (I in Fig. ^ drops to zero at 
a certain critical value of the confinement strength (neglecting Zeeman splitting). Both the 
value for the critical confinement strength and the overall shape of the curve agree well with 
the results of Ref. ||. 

An important consequence of exchange on the distribution of electrons for experiments 
is the increased separation between the compressible regions. This causes a strong decrease 
in inter-edge channel scattering because of the reduced overlap of the edge-channel wave 
functions, which leads to an increase in the spin-flip equilibration length of the edge chan- 
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nels. This in contrast to the assumption of an edge channel separation of the order of 
one magnetic length in Ref. Using a typical edge channel separation of a few magnetic 
lengths in accordance with the present results, Khaetskii's theory0 would give a much longer 
equilibration length. 

The exchange enhanced channel separation has also consequences for the two-terminal 
magnetoconductance of a narrow channel or point contact. In Ref. ^, Chklovskii et al. 
discuss the conductance for these systems in the framework of their electrostatic descrip- 
tion of the channels. They obtain conductance quantization but, due to the small width 
of the incompressible regions, the calculated plateau widths are much smaller than those 
experimentally observed. They attribute this discrepancy between theory and experiment to 
the presence of disorder. However, we propose as an alternative explanation the electronic 
exchange interaction which we find to strongly enhance the width of the insulating strips. 



B. Accuracy 

In this subsection we study the accuracy of our theory and how our present results are 
altered if we take into account realistic features which do not change the qualitative results 
but are important to describe experimental situations. 

First we investigate the validity of Eq. ([TtD for high magnetic fields where mixing with 
higher Landau levels was disregarded. To this end we use a parabolic confinement \m*ujQX^ 
and a density distribution of the form: n{x) = — W/2)9{x + W/2), where W is the 

width of the wire. The 2DEG in the presence of the uniform magnetic field and parabolic 
confinement potential is our unperturbed system and we treat the Hartree potential as a 
perturbation. The energy levels of the unperturbed system are Enk = {n+\)hn + h'kl/2M, 
where Vt = ^u'^ + Uq is the renormalized cyclotron frequency and M is the renormalized 
electron mass: M = m*{l + ujI/uJq). The wave functions are the bulk Landau level wave 
functions with effective magnetic length I = Ib\JojI/{uoI + c<Jq).0 According to first order 
perturbation theory, the ground state wave function, of the electron at the edge, with 
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quantum number X = W/2, is approximately: 

m = m) + —E^ — ^ — l^i), (20) 

where we have used the Hartree potential as the perturbation. A straightforward calculation 
of the coefficient of {ipi) (which of course depends on W) shows that for 20 /b < < 40 /s, 
a magnetic field of 5 T and a typical value of hujQ = 3 meV for the parabolic confinement, 
the mixing with the first Landau level wave function is approximately 10% and that this 
value decreases as 1/\^Bq with increasing magnetic field. We have to conclude that the dip 
in the electrostatic, i.e. Hartree and confinement, potential at the edges which is responsible 
for the spin-polarizing transition and ultimately for the compressible edge state, cannot be 
compensated by mixing with higher Landau levels, i.e. that our approximation is allowed 
in the high magnetic field limit considered here. By performing a similar calculation it is 
possible to show that the mixing due the exchange-correlation potential used by Ferconi and 
Vignale0 is only a few percent. 

Another approximation which gives rise to quantitative deviations from our theory is the 
assumption that the electron gas is strictly two dimensional. In reality the electron wave 
functions also extend in the 2;-direction and a form factor F{q) in the Coulomb interaction, 
V{q) = e^F(g)/2e|g|, takes this effect into account. Calculations of edge channel splitting 
including a form factor in the Hartree-Fock approximation have already been calculated by 
Rijkels and Bauer .0 By comparison with Ref. ^we estimate the values of a to decrease by 
about 0.004 if a form factor would be included into our calculations. This shifts the range 
of hujQ in which both the electrostatic and the Hartree-Fock regime are unreliable to lower 
energies, i.e. to 3.2 meV-5.1 meV. 

IV. APPLICATIONS 

The present theory can be easily generalized to a variety of systems. In the next subsec- 
tions we consider two special cases, which are experimentally relevant. 
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A. Back gate 



As a first application we have carried out calculations for a system like the one considered 
in Sec. |11 A|, but now in the presence of an infinite plane of constant potential at a distance d 



from the 2DEG, which may represent a real back gate or a plane of not fully depleted donors 
parallel to the 2DEG. This gate can be easily included into the calculations by adding the 
potential of a mirror charge distribution of opposite sign located at a distance 2d from the 
electron gas: 

Vf,g{x) = —^ r dx n{x) In - x'f + {2df . (21) 

In Fig. ^ we have plotted the solutions that minimize the total energy of the system including 
the back gate as a function of the distance d to the wire for a given number of electrons and 
a fixed confinement potential. 

It is clear from the figure that moving the back gate closer to the 2DEG has the effect of 
reducing the width of the compressible regions. This is to be expected because the back gate 
screens the long-range Coulomb interaction in the electron gas, thus effectively reducing the 
importance of the direct interaction relative to the exchange interaction and thereby forcing 
the solution for small distances d into the Hartree-Fock regime. 

It should be noted that the description of the back gate given here does not apply for 
a quantum wire defined by a depletion gate. In that case the back gate would also screen 
the charge on the depletion gates, thus increasingly altering the confinement potential as 
the back gate is moved closer to the 2DEG. In that case a fully self-consistent calculation 
is necessary. Furthermore, for our results to be correct, the distance from the gate to the 
2DEG should not be too small since in that case also the short-range exchange interaction 
would be influenced by the gate. 
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B. Ground state of a quantum dot 



A system that is very similar to the strip previously considered is a quantum dot in a 
high magnetic field, especially when the number of electrons in the dot is large. There are 
two main differences. In the first place we use the symmetric gauge to describe the dot. In 
this gauge the wave functions in the lowest Landau level, labeled by the quantum number 
of angular momentum, are (f)m{z) = where z = The second difference is 

that because the system is finite the filling factor u is no longer labeled by the continuous 
variable X but by the discrete quantum number m. As a consequence the radial symmetric 
density is now given by: 

where u!^ is the filling factor for a given angular momentum m and spin direction a. We 
impose as boundary conditions that < z/^ < 1 and that the total charge in one spin level 
J2m is integer. Due to the circular symmetry of the dot the Hartree potential is in 
this case given by: 

P rliv roc t'tjIt'i 

Vn{r) = --^ d<p dr' / ^ (23) 



An analysis analogous to the one presented in Sec. [111B| shows that for high magnetic 
fields (approximately 5 T) the mixing with higher Landau levels due to the Hartree po- 



tential and the exchange- correlation potentials used in Ref. |TI| is only a few percent. As a 
consequence, fractional occupation numbers are essential, which cannot be described by the 



Kohn-Sham scheme of Ref. 11 



We have calculated the ground state of a dot consisting of 40 electrons using a parabolic 
confinement potential Vc{r) = a ^^^^^ with varying strength a in a magnetic field of 5.0 T. 
To find the ground state for a given magnetic field we vary the distribution of electrons among 
the occupied spin levels while keeping the total number of electrons constant. The results 
are plotted in Fig. M. The ground states are similar to those of the strip except for the fact 
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that the middle region is wider. This is due to the fact that the wave functions are centered 
around the radii = V2m Is, where m is the angular momentum index, and not around 
equidistant points like in the strip. 

In Ref. a selfconsistent calculation was performed for the density profile of a quantum 
dot containing 39 electrons using a modified Hartree form for the electrostatic electron- 
electron interaction: 

Ke(r, r') = --^ I ^ = - ^ I , (24) 

thus taking into account the finite z-extent {6z) of the wave functions and a back gate at 
distance d. Again we want to extend these calculations to include exchange but we must 
keep in mind that a finite {6z) also reduces the exchange interaction by a form factor: 

F{{6z)) = exp (^(5z)^)Erfc (^^^ , (25) 

where {6z) is given in units of the magnetic length and Erfc(x) = 1 — Erf(x) denotes the 
complementary error function. A calculation using interaction Eq. ( p^ and the local ex- 
change contribution reduced by the form factor Eq. (^5]) gives a ground state which differs 
from that of McEuen et al. in the sense that the widths of the incompressible (compressible) 
regions are larger (smaller) in our solution. 

We have calculated the addition spectrum of dots consisting of 37 and 38 electrons in 
order to determine the effect of exchange on the inter-level transitions that are responsible 
for the oscillatory behavior of the addition energy. Here we take into account the fact that 
the magnetic length is renormalised by the magnetic field, Eq. (0), in order to compare 
the energies at different magnetic fields. The addition energy /S.E{N) needed for adding an 
extra electron to a dot of electrons is: 

^E{N) = E{N + l)-E{N), (26) 

where E{N) is the total energy given by Eq. (|^). The results are depicted in Fig. |^. It 
should be mentioned that our results are very close to the Hartree-Fock solutions. This is 
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due to the fact that the combination of the back gate and exchange force the solution into 
the integer fiUing regime. 

The inset of Fig. |^ shows the energy dependence of different dot configurations with 
magnetic field for a dot containing 39 electrons. The ground state is formed by the config- 
uration with lowest energy. The numbers in the inset represent the number of electrons in 
the majority {i.e. lowest) spin level for the indicated curve. As the magnetic field increases, 
the degeneracy of the spin levels is increased and at certain magnetic fields an electron can 
jump from the upper to the lower level, causing the kinks in the ground state energy. When 
all electrons are in the lowest spin level no transitions can occur and the curves for both the 
total and addition energy are smooth. 

Recent calculationd^ and measurementsil indicate that at higher magnetic fields than the 
ones shown in Fig. ^ an edge reconstruction occurs, resulting in an extra kink in the curve 
for the addition energy in the regime where only one spin level is occupied. Performing a 
calculation with filling factors that are restricted to integer values, we are able to reproduce 
this feature. However, for the small confinement strength a at which this transition occurs, 
the fractional-filling solution is not so well converged. In spite of the fact that we find 
evidence for the edge reconstruction, the numerical accuracy of our solution does not allow 
decisive conlusions. 



In comparison with the results obtained in Ref. the oscillations found in our calcu- 
lations have a larger period and amplitude by a factor of approximately 1.5. The increase 
in amplitude of the oscillations can be qualitatively understood in terms of the capacitance 
model for the island proposed by Evans et a/.ii and is due to the increase of the width of 
the incompressible strip in the picture including exchange. This results in a decrease of the 
capacitance between the inner and outer compresible regions where the extra electron can 
be added and hence in an increase in the peak-to-peak amplitude. The larger period can 
be explained in a similar fashion: Exchange effects reduce the widths of the compressible 
strips and thus the capacitances Ci and C2 between the respective compressible strips and 
the back gate which results in an increased period of the oscillations compared to the elec- 



trostatic treatment of McEuen et al. Compared with their experimental data, the period 
of the oscillations found in our calculations is about 0.3 T too large whereas their calcu- 
lations give a period which is 0.25 T too small. However, these values depend sensitively 
on the calculated ground state density since the total energy, with a typical value around 
1 eV, has to be calculated with an absolute accuracy of at least 0.1 meV in order to resolve 
the oscillations and therefore small numerical deviations from the true ground state would 
drastically influence the addition energies. 



V. SUMMARY AND CONCLUDING REMARKS 

We have used density-functional theory in strong magnetic fields to investigate the ef- 
fects of the electron-electron interaction on edge states in the integer quantum Hall regime. 
We have included exchange in the local density approximation. In this approximation 
smooth density distributions corresponding to fractional filling at the edges can be treated, 
which is beyond a Hartree-Fock treatment of edge channels.^ To describe fractional occupa- 
tion could otherwise only be achieved by exact diagonalization or configuration interaction 
calculations.! We have found that the width of the edge channels is strongly reduced due 
to exchange effects and that the width of the incompressible strips is of the order of several 
magnetic lengths, which should have a strong influence on the spin-flip equilibration length 
and the two-terminal conductance of a point contact. We predict a range of confinement 
potentials huJo=3.2 meV-5.1 meV in which the entire regime from electrostatic to Hartree- 
Fock is covered. Furthermore we have calculated ground state properties of the system in 
the presence of a back gate parallel to the quantum wire. We established an increasingly 
important role for the exchange interaction as the back gate is closed in on the 2DEG. As 
a second application we have calculated the ground state of a quantum dot in a strong 
magnetic field. We found that the results are similar to those found for the quantum wire. 
We compared addition energies for the dot with the self consistent calculations of McEuen et 
a/.0 and found that the oscillations in the addition energies have larger amplitudes and 
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periods in our calculations. 
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FIGURES 

FIG. 1. Schematic picture of the system, consisting of a strip of 2DEG in a uniform magnetic 
field. The wire is confined by a parabohc potential in the x-direction. For details sec text. 

FIG. 2. Electronic ground state filling factors, (a) ignoring exchange and including exchange 
in (b) the Hartree-Fock approximation and (c) the local density approximation. The filling fac- 
tors are plotted for confinement strengths a = 0.047 (solid line), 0.041 (dashed line) and 0.035 
(dashed-dotted line). Calculated for a magnetic field of 7.2 T. 

FIG. 3. One-particle potentials and density distribution for the ground state of Fig. 2(c) with 
a = 0.035 for (a) majority-spin electrons and (b) minority-spin electrons. Plotted are the electro- 
static potential consisting of the confinement and the Hartree potential (dashed line) , the exchange 
potential (dashed-dotted line) and the total potential excluding the Zeeman term (solid line). In 
the regions where the exchange potential is zero the total and Hartree potential coincide. The 
potentials are offset for clarity and are given in units of ^^^^^ . 

FIG. 4. Width of innermost compressible region (C) and outermost incompressible region (I) 
(see inset), plotted as a function of the confinement strength a. The solid lines are drawn as a 
guide for the eye. 

FIG. 5. Occupation of the system in the presence of the back gate at a distance d = oo 
(dashed-dotted line), d = 201b (dashed line) and d = 101b- The confinement strength a = 0.035 
and B = 7.2 T. (solid line). 

FIG. 6. Occupation numbers of a quantum dot consisting of 40 electrons for confinement 
strengths a = 0.09 (solid line), 0.08 (dashed line) and 0.07 (dashed-dotted line). The discrete 
filling factors are positioned at the radii and have been connected by lines for clarity. 
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FIG. 7. Addition energy for quantum dots of 37 and 38 electrons as a function of magnetic 
field. The calculation was performed for a confinement strength Hloq = 1.6 meV, d = 100 nm and 
(dz) = 10 nm. The inset shows the ground state energies of the different possible configurations 
of the dot as a function of magnetic field. The numbers represent the number of electrons in the 
lowest spin level. Only odd numbers are indicated. 
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